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Abstract 

In 1968 we published the solution of the ground state energy and 
wave function of the one-dimensional Hubbard model, and we also 
showed that there is no Mott transition in this model. Details of the 
analysis have never been published, however. As the Hubbard model 
has become increasingly important in condensed matter physics, re- 
lating to topics such as the theory of high-T c superconductivity, it is 
appropriate to revisit the one-dimensional model and to recall here 
some details of the solution. 
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1 Introduction 

In a previous paper QTJ we reported the solution of the one- dimensional (ID) 
Hubbard model, showing the absence of the Mott transition in its ground 
state, but the letter format of the paper did not permit the presentation 
of all the details of the analysis. Over the years the Hubbard model [§, || 
has become more important, for it plays an essential role in several topics 
in condensed matter physics, including ID conductors and high-T c supercon- 
ductivity. It also plays a role in the chemistry of aromatic compounds (e.g., 
Benzene @, j5|). Several books || [7|, || || now exist in which the ID Hubbard 
model is analyzed, and numerous papers have been written on properties of 
the model [Hj. Almost invariably these publications are based upon results of 



|[J], including the absence of a Mott transition, but without derivation. While 
other rigorous results on higher dimensional Hubbard models exist, and some 



of these are reviewed in flll], |12[, the ID model stands as the only Hubbard 
model whose ground state can be found exactly. It has been brought to our 
attention that it would be useful to students and researchers if some details 
of the solution could be made available. Here, taking the opportunity of the 
symposium, StatPhys- Taiwan 2002, which takes place in the year when both 
of the authors turn 70, we revisit the ID Hubbard model and present some 
details of the 34-year old solution. 

While our paper |l| contained significant results about the excitation 
spectrum, it was mainly concerned with the integral equations for the ground 
state and we concentrate on those equations here. The new, unpublished 
results are contained in Sects. 5, 6, 7. 

Consider a crystal of N a lattice sites with a total of N itinerant electrons 
hopping between the Wannier states of neighboring lattice sites, and that 
each site is capable of accommodating two electrons of opposite spin, with 
an interaction energy U > 0, which mimics a screened Coulomb repulsion 
among electrons. The Hubbard model is described by the Hamiltonian 

H = T ^ Yl c ^ c i ° + U Yl Ui T n n , (!) 

<ij> a i 

where c\ a and q„ are, respectively, the creation and annihilation operators 
for an electron of spin a in the Wannier state at the z-th lattice site and 
n ia = c\ a Ci a is the occupation number operator. The summation < ij > 
is over nearest neighbors, and one often considers (as we do here) periodic 



3 



boundary conditions, which means that < ij > includes a term coupling 
opposite edges of the lattice. We are interested in the ground state solution 
of the Schrodinger equation TC = E . 

For bipartite lattices (i.e., lattices in which the set of sites can be divided 
into two subsets, A and B, such that there is no hopping between A sites or 
between B sites), such as the one-dimensional chain, the unitary transfor- 
mation V^TLV leaves TC unchanged except for the replacement of T by — T. 
Here V = exp [in J2 i£A (rii ^ + rii , with A being one of the two sublattices. 
Without loss of generality we can, therefore, take T = —1. In any event, 
bipartite or not, we can renormalize U by redefining U to be U/\T\. Hence- 
forth, the value of T in (|l|) is —1 and U is positive and fixed. The dependence 
of the Hamiltonian and the energy on U will not be noted explicitly. 

The commutation relations 



rii t, TC 







imply that the numbers of down-spin electrons M and up-spin electrons M' 
are good quantum numbers. Therefore we characterize the eigenstates by M 
and M', and write the Schrodinger equation as 

H\M,M') = E(M,M')\M,M'). (2) 

Naturally, for any fixed choice of M, M' there will generally be many solutions 
to (0), so that \M,M') and E(M,M') denote only generic eigenvectors and 
eigenvalues. Furthermore, by considering particles as holes, and vice versa, 
namely, introducing fermion operators 

d- = A d ] = c- 

and the relation rii a — 1 — d] a di a , we obtain the identity 

E(M, M') = ~(N a - N)U + E(N a — M,N a — M') (3) 

where 

N = M + M' 

is the total number of electrons. Since N > N a if, and only if, (M — A^ a ) + 
(M' — A^ a ) < N a , we can restrict our considerations to 



N<N a , 
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namely, the case of at most a "half-filled band". In addition, owing to the 
spin-up and spin-down symmetry, we need only consider 

M < M'. 



2 The one-dimensional model 

We now consider the one-dimensional model, and write |M, M') as a linear 
combination of states with electrons at specific sites. Number the lattice sites 
by 1, 2, • - • , N a and, since we want to use periodic boundary conditions, we 
require N a to be an even integer in order to retain the bipartite structure. 
For later use it is convenient also to require that N a = 2 x (odd integer) in 
order to be able to have M = M' = N a /2 with M odd. For the ID model 
the sum in ([]]) over < ij > is really a sum over 1 < % < N a ,j = i + 1, plus 
1 < j < N a , % = j + 1 with iV + 1 = 1. 

Let | Xi, ■ ■ ■ ,xn) denote the state in which the down-spin electrons 
are located at sites x%, x 2 , ■ ■ ■ ,Xm an d the up-spin electrons are at sites 
x m+i, " " " j^jv- The eigenstate is now written as 

\M,M')= Yl f(xi, ■■■ ,x N )\x u ■■■ } x N ) (4) 

l<Xi<N a 

where the summation is over all x\, • • • , xn from 1 to N a , and f(x\, • • • , xn) 
is the amplitude of the state | x±, ■ ■ ■ ,Xjy). 

It is convenient to denote the iV-tuple Xi,x 2 , ...,Xn simply by X. 

By substituting (|) into the Schrodinger equation (|j), we obtain (recall 
T = -1) 

N 

^ [/(Xi, • • • ,Xj + 1, • • • ,X N ) +/(gi, • • • ,Xj - 1, • • • ,X N ) 

1 

u\^2S(x i -x j ) f(xi, ■ ■ ■ ,x N ) = E f(xx, ■ ■ • ,x N ), (5) 



where 5 is the Kronecker delta function. We must solve (^) for / and E, with 
the understanding that site is the same as site N a and site N a + 1 is the 
same as site 1 (the periodic boundary condition). Equation @ is the 'first 
quantized' version of the Schrodinger equation (fj). It must be satisfied for 



all 1 < Xi < N a , with 1 < i < N. 
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As electrons are governed by Fermi-Dirac statistics, we require that f(X) 
be antisymmetric in its first M and last M' variables separately. This anti- 
symmetry also ensures that / = if any two x's in the same set are equal, 
which implies that the only delta-function term in (|5]) that are relevant are 
the ones with i < M and j > M. This is consistent with the definition of 7i 
in (0), in which the only interaction is between up- and down-spin electrons. 

The antisymmetry allows us to reinterpret fl5|) in the following alternative 
way. Define the region R to be the following subset of all possible values of 
X (note the < signs): 

/ 1 < acx < x 2 < ■ ■ ■ < x M < N a \ | , 

V 1 < X M+1 < XM+2 < ■■■ <X N <N a J J • 

In R any of the first M Xj's can be equal to any of the last M', with an 
interaction energy nil, where n is the number of overlaps of the first set with 
the second. 

The antisymmetry of / tells us that / is completely determined by its 
values in the subset R, together with the requirement that / = if any 
two s's in the same set are equal (e.g., X\ = x 2 ). Therefore, it suffices to 
satisfy the Schrodinger equation ([]) when X on the right side of (|5|) is only 
in R, together with the additional fact that we set / = on the left side 
of (|D if Xi ± 1 takes us out of R, e.g., if x± + 1 = x 2 - (Warning: With 
this interpretation, equation (P) then becomes a self-contained equation in 
R alone and one should not ask it to be valid if X R.) 

There is one annoying point about restricting attention to R in When 
X\ — 1 the left side of (||) asks for the value of / for x\ — N a , which takes us 
outside R. Using the antisymmetry we conclude that 

/(JV a , x 2 , x N ) = (-l) M ~ l f(x 2 , ...,x M , N a , x M +i, x N ), (7) 

with similar relations holding for xm = 1, %m+i = N a or xjy = 1. Equation 
(0) and its three analogues reflect the "periodic boundary conditions" and, 
with its use, ([5[) becomes a self-contained equation on R alone. 

We now come to the main reason for introducing R. Let us assume that 
M and M' are both odd integers. Then (— 1) M_1 = (— 1) M _1 = 1 and we 
claim that: For all U , the ground state of our Hamiltonian satisfies 

1. There is only one ground state and 

2. f(X) is a strictly positive function in R. 
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To prove 2. we think of (|^) as an equation in R, as explained before. 
We note that all the off-diagonal terms in the Hamiltonian (thought of as a 
matrix 7i from £ 2 (R) to £ 2 (R)) are non-positive (this is where we use the fact 
that (-l)^- 1 = (-l)^'- 1 = 1). If E Q is the ground state energy and if f(X) 
is a ground state eigenfunction (in R), which can be assumed to be real, then, 
by the variational principle, the function g(X) = \ f(X)\ (in R) has an energy 
at least as low as that of /, i.e., (g\ g) = (f \ f) and (g\TC\g) < (f\TC\f) since 
\f(X)\H(X,Y) \f(Y)\ < f(X)H(X,Y)f(Y) for every X,Y. Hence g must 
be a ground state as well (since it cannot have a lower energy than Eq, by 
the definition of Eq). Therefore, g(X) must satisfy (§) with the same Eq. 
Moreover, we see from ([5]) that g(X) is strictly positive for every X G R 
(because if g(Y) = for some Y G R then g(Z) = for every Z that differs 
from Y by one 'hop'; tracing this backward, g(X) = for every X G R). 

Returning now to /, let us assume the contrary of 2., namely, f(X) > 
for some X G R, and f(Y) < for some Y G R. We observe that since 
h = g — f must also be a ground state (because sums of ground states are 
ground states, although not necessarily normalized), we have a ground state 
(namely h) that is non-negative and non-zero, but not strictly positive; this 
contradicts the fact, which we have just proved, that every no n- negative 
ground state must be strictly positive. Thus, 2. is proved. 

A similar argument proves 1. If / and /' are two linearly independent 
ground states then the state given by k(X) = f(X) + cf'(X) is also a ground 
state and, for suitable c, k(X) = for some X G R, but k cannot be 
identically zero. Then \k\ will be a non- negative ground state that is not 
strictly positive, and this contradicts statement 2. 

The uniqueness statement 1. is important for the following reason. Sup- 
pose that we know the ground state for some particular value of U (e.g., 
U = oo) and suppose we have a [/-dependent solution to (§) in some inter- 
val of U values (e.g., (0, oo)) with an energy E(U) such that: (a) E{oo) is 
the known ground state energy and (b) E{U) is continuous on the interval. 
Then E(U) is necessarily the ground state energy in that interval. If not, the 
curve E(U) would have to cross the ground state curve (which is always con- 
tinuous), at which point there would be a degeneracy - which is impossible 
according to 1. 

Items 1. and 2. can be used in two main applications. The first is the 
proof of the fact that when M and M' are odd the ground state belongs 
to total spin S equal to \M — M'\/2 and not to some higher S value. The 
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proof is the same as in [13]. In [13] this property was shown to hold for all 
values of M and M', but for an open chain instead of a closed chain. In the 
thermodynamic limit this distinction is not important. 

The second main application of these items 1. and 2. is a proof that the 
state we construct below using the Bethe Ansatz really is the ground state. 
This possibility is addressed at the end of Sect. |3] where we outline a strategy 
for such a proof. Unfortunately, we are unable to carry it out and we leave 
it as an open problem. 

We also mention a theorem [14], which states that the ground state is 
unique for M = M' = N a /2 and N a = even (the half-filled band). There is 
no requirement for M = M' to be odd. 



3 The Bethe Ansatz 



The Bethe Ansatz was invented |L5] to solve the Heisenberg spin model, which 
is essentially a model of lattice bosons. The boson gas in the continuum 
with a positive delta function interaction and with positive density in the 



thermodynamic limit was first treated in ||16|| . McGuire []T7| was the first 
to realize that the method could be extended to continuum fermions with a 
delta function interaction for M — 1. (The case M = is trivial.) The first 
real mathematical difficulty comes with M = 2 and this was finally solved in 
T8|. The solution was inelegant and not transparent, but was a precursor to 



the full solution for general M by Gaudin [|T^] and Yang |2C 



We now forget about the region R and focus, instead, on the fundamental 
regions (note the < signs) 

R Q = {X : 1 < x Q1 < ■ ■ ■ < x QN < N a } (8) 

Here Q = {Ql, Q2, • • • , QN} is the permutation that maps the ordered set 
{1,2, ...,N} into {Ql, Q2, QN}. There are iV! permutations and corre- 
sponding regions Rq. The union of these regions is the full configuration 
space. These regions are disjoint except for their boundaries (i.e., points 
where x Qi = x Q ( i+1 )). 

Let k\ < k 2 < ... < A; at be a set of unequal, ordered and real numbers 
in the interval — tc < k < it, and let [Q, P] be a set of iV! x iV! coefficients 
indexed by a pair of permutations Q, P, all yet to be determined. 
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When X 6 Rq we write the function f(X) as (the Bethe Ansatz) 

/(X) = fq(xt, ■ ■ ■ ,x N ) = [Q,P] exp[i(k P1 x Q1 H h k PN x QN )] . (9) 

p 

In order for (|S|) to represent a function on the whole configuration space 
it is essential that the definitions @ agree on the intersections of different 
Rq's. This will impose conditions on the [Q, P]'s. 

Choose some integer 1 < z < AT and let j = i + 1. Let P, P' be two 
permutations such that Pi = P'j and Pj = P'i, but otherwise Pm = P'm 
for m 7^ i, j. Similarly, let Q,Q' be a pair with the same property (for this 
same choice of i) but otherwise P, P' and Q, Q' are unrelated. 

The common boundary between Rq and Rq> is the set in which xq-i = xqj. 
In order to have /q = /q/ on this boundary it is sufficient to require that 

[Q 1 P} + [Q 1 P'} = [Q',P} + [Q' 1 P'] . (10) 

The reason that this suffices is that on this boundary we have xqi = xqj and 
kpi + k P j = kpq + kp/j. Thus, fllCf ) expresses the fact that the exponential 
factor exp[i(k Pi XQi + k P jXQj)} = exp[i(kp'iXQi + k P >jXQj)} is the same for Q 
and Q', and for all values of the other x m 's. 

Next we substitute the Ansatz (Q) into (|5|). If \xi — Xj\ > 1 for all i,j 
then, clearly, we have 

N 

E = E{M,M') = -2 cos k i- ( n ) 

i=i 

We next choose the coefficients [Q, P] to make (0) hold generally — even 
if it is not possible to have |xj — xA > 1 for all i,j when the number of 
electrons exceeds N a /2. The requirement that (|TT]) holds will impose further 
conditions on [Q,P] similar to flip]) . 

Sufficient conditions are obtained by setting XQi = XQj on the right side 
of (H) and requiring the exponential factors with xq^ and XQj alone to satisfy 
d^). In other words, we require that 

[Q, P] e~ ikp J + [Q, P'] e- ik *"i + [Q,P) e +lkp > + [Q, P'] e +lkp '> 
= [£', P] e _i *" + [Q 1 , P'} e-' lkp '< + [Q', P] e +ife ^' + [Q', P'] e +ik *"i 

+U([Q,P] + [Q,P'\) . (12) 
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If we combine ([12]) with (0) and recall that kpj = kp'i, etc., we obtain 
= • I ~ %U i 2 rr779 ^ p/ 1 

sm fcpj — sin kpj + lU / 2 

+ sinfcp, -sin fc Pj - ^ 
sin kpi — sin fcpj + iU/2 

It would seem that we have to solve both (pT3|) and ([K]) for the (AH) 2 
coefficients [Q,-P], and for each 1 < i < N — 1. Nevertheless, ( |13|) alone is 
sufficient because it implies (|T0|). To see this, add ([T3"D, as given, to ( pT3|) with 
[Q', -P] on the left side. Since Q" = Q, the result is (TTOf) . Our goal, then, 
is to solve (|I3|) for the coefficients [Q, P] such that the amplitude / has the 
required symmetry. 

These equations have been solved in |0J and [^TJ, as we stated before, 
and we shall not repeat the derivation here. In these papers the function 
sin k appearing in ([13]) is replaced by k, which reflects the fact that JOJ and 



20| deal with the continuum and we are working on a lattice. This makes no 
difference as far as the algebra leading to equations fll4j) is concerned, but it 
makes a big difference for constructing a proof that these equations have a 
solution (the reason being that the sine function is not one-to-one). 

The algebraic analysis in [05] and pO] leads to the following set of iV + M 
equations for the N ordered, real, unequal fe's. (Recall that M < M'.) They 
involve an additional set of M ordered, unequal real numbers Ai < A2 < 
• • ■ < A M . 

m rr ^sin kj — iAg — U/4 , , 

^ 1 sm kj — iKp + U /4 

n V isin kjj-iAg -U/4 = _ y4 -zAg + iA a + ?7/2 = 
i sin fc, - zA a + U/4 J- A -iAfl + zA a - C//2 ' " ' ' ' " ' 

J=l J /5=1 K 

We remark that an explicit expression for the wave function f(X) has been 
given by Woynarovich j^TJ, part 1, Eqs. 2.5-2.9]. 

These equations can be cast in a more transparent form (in which we now 
really make use of the fact that the fc's and A's are ordered) by defining 

0(p) = -2 tan- 1 (^p) , -71 < 9 < tt 
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Then, taking the logarithm of (|14|), we obtain two sets of equations 



N a k j = 2nI j + ^29(2smk j -2Ap), j = 1, 2, • • • ,N (15) 

N M 

^0(2sin ^-2A a ) =27rJ Q -^0(A a -A /3 ), a = l,2,--- ,M (16) 

j=l (3=1 

where Ij is an integer (half-odd integer) if M is even (odd), while J a is an 
integer (half-odd integer) if M' is odd (even). 

It is noteworthy that in the U — > oo limit the two sets of equations 
essentially decouple. The A's are proportional to U in this limit, but the sum 
in (|T5D becomes independent of j. In particular, when the A's are balanced 
(i.e., for every A there is a —A) as in our case, then this sum equals zero. 

From (|15|) and ( Jl6| ) we have the identity 



N n N M 

Z7T 



i ~ a -i -i 



For the ground state, with = 2x (odd integer) and M = N/2 = odd, we 
make the choice of the Ij and J a that agrees with the correct values in the 
case U = oo, namely 

Ij = j — (N + l)/2 J a = a - (M + l)/2 (18) 

We are not able to prove the existence of solutions to fll5| ) and ( |16|) that 
are real and increasing in the index j and a. In the next section, however, 
we show that the — >• oo limit of ([15]) and fllfiD has a solution, and in Sect. 
^| we obtain the solution explicitly for N/N a = 1. This leaves little doubt 
that ( |l~5|) and flTB]) can be solved as well, at least for large N. 

Assuming that M = M' = N/2 is odd, the solution is presumably unique 
with the given values of Ij and J a and belongs to total spin S = 0. 

Assuming that the solution exists, we would still need a few more facts 
(which we have not proved) in order to prove that the Bethe Ansatz gives 
the ground state. 

(a.) Prove that the wave function @ is not identically zero. 

(b.) Prove that the wave function (||) is a continuous function of U. 
From the uniqueness of the ground state proved in Sect. |2], and the fact 
that the solution (0) coincides with the exact solution for U — oo (in which 



11 



case Jq{x) is a Slater determinant of plane waves with wavenumbers kj = 
2iiIj/N a ), (a.) and (b.) now establish that the wave function (§) must be 
the ground state for all U. 

Remark: Assuming that the Bethe Ansatz gives the ground state for a 
given M < M' then, as remarked at the end of Sect. 2 (and assuming M and 
M' to be odd) the value of the total spin in this state is 5* = (M' — M)/2. 
Thus, the solution to the Bethe Ansatz we have been looking at is a highest 
weight state of SU(2), i.e., a state annihilated by spin raising operators. 



4 The ground state 

For the ground state Ij = I(kj) and J a = J(A a ) are consecutive integers or 
half-odd integers centered around the origin. As stated in Sect. |3], each kj 
lies in [— n, it] (since kj — *> kj + 2im defines the same wave function). In the 
limit of N a , N, M, M' — > oo with their ratios kept fixed, the real numbers k 
and A are distributed between — Q and Q < ir and — B and B < oo for some 
< Q < 7i and < B < oo. In a small interval dk the number of k values, 
and hence the number of j values in (|15j), is N a p(k)dk, where p is a density 
function to be determined. Likewise, in a small interval dA the number of A 
values and a values in ( |TED is N a a(A)dA. An alternative point of view is to 
think of I(k) as a function of the variable k. Then I(k + dk) — I(k) counts 
the number of k values between k and k + dk so we have dJ(fc) /dk = N a p{k) . 
A similar remark holds for J (A). 

The density functions p(k) and cx(A) satisfy the obvious normalization 

p(k)dk = N/N a , [ a{A)dA = M/N a . (19) 

J J-B 

By subtracting (|T5| ) with j from (|15|) with j + N a p(k)dk, and taking the 
limit N a — ► oo we obtain ([2t]) below. Likewise, subtracting ( JIB| ) with a from 



(16) with a + N a a(A)dA, and taking the limit N a — >• oo we obtain (pT|). 
An alternative point of view is to take the derivatives of (|T5|) and flTB| ) with 
respect to kj and A a respectively, set dl /dk = N a p(k), dJ/dA = N a a(A), 
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and take the N„ — » oo limit. 



1 = 2tt p( k) + 2 cos k [ dA a (A) 6' (2 sin k - 2A) 

J-B 

rQ 

-2 / dkp(k)9' (2 sin k-2A) 

J-Q 



B 

2vra(A)- / dA' a(A') 9' (A - A') 



or, equivalently, 



where 



A(A - A') 
A^ 2 (A - A') 



--0'(2A-2A') = — 

— e'(A-A') = — 

2tt v y 2tt 



817 



L[/ 2 + 16(A - A') 2 
AU 



U 2 + 4(A - A') 2 

/oo 
K(A — x)K(x — A') dx. 
-oo 



(20) 




(21) 



(22) 



(23) 



Note that .fT 2 is the square of K in the sense of operator products. Note also 
that (p2|) and (|23f ) are to be satisfied only for \k\ < Q and |A| < B. Outside 
these intervals p and a are not uniquely defined, but we can and will define 
them by the right sides of (|2|) and (p3|) . 

The following Fourier transforms will be used in later discussions: 



/oo /*oo 
e^ A A-(A) dA = e ~ u ^\ / e ia,A iT 2 (A) dA = e~ u ^l 2 . (24) 
•oo J — oo 

The ground state energy (|TTD now reads 

rQ 

E(M, M') = —2N a / p(k) cos fc dk. (25) 
</-<2 



13 



where p(k) is to be determined together with cx(A) from the coupled integral 
equations (p2|) and (|23|) subject to the normalizations ("|1QD . 

5 Analysis of the integral equations 



In this section we shall prove that equations (j22|) and (^3|) have unique solu- 
tions for each given < Q < tt and < B < oo and that the solutions are 
positive and have certain monotonicity properties. These properties guaran- 
tee that the normalization conditions ( |T9| ) uniquely determine values of Q, B 
for each given value of N when M = M' = N/2 (in this case we have B = oo). 
However, we have not proved uniqueness of Q, B when M ^ M' (although 
we believe there is uniqueness). But this does not matter for the absolute 
ground state since, as remarked earlier, the ground state has S = (in the 
thermodynamic limit) and so we are allowed to take S z = 0. For M ^ M', 
we have remarked earlier that the solution probably has S = \M' — M\/2 
and is the ground state for S = \M' — M\/2. 

An important first step is to overcome the annoying fact (which is relevant 
for Q > 7r/2) that sin k is not a monotonic function of k in [— tt, n]. To do this 
we note that (cos k)K (sm k — A) is an odd function of k — n/2 (for each A) 
and hence p(k) — 1/2-k also has this property. On the other hand, i^(sin k, A) 
appearing in (p3|) is an even function of k — tt/2. As a result p(k) appearing 
in the first term on the right side of (B3f) can be replaced by 1/2tt in the 



intervals Q' < k < Q and — Q < k < —Q', where Q' = n — Q. Thus, when 



Q > tt/2, we can rewrite the [Q', Q] portion of the first integral in (p3|) as 
[ Q K(smk-A)p(k)dk= ^ K(sin k - A)— dk = — I ' K(sink — A) 

Jq> Jq> 27T 2TT Jq, 



dk 



A similar thing can be done for the [-Q, —Q'] portion and for the corre- 
sponding portions of (O). 

The integrals over k now extend at most over the interval [—tt/2, tt/2], in 
which sin k is monotonic. 

We are now in a position to change variables as follows. For — 1 < x < 1 

let 

t( x ) = -*-(l - x 2 )~ 1/2 , f(x) = (1 - ^-^(sin" 1 x) . (26) 
2tt 
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In case Q < ir/2, p(sin~ 1 x) is denned only for sinx < Q, but we shall soon 
see (after ( ^8|) how to extend the definition of / in this case. We define the 
step functions for all real x by 

B(x) = 1, \x\ < B ; = 0, otherwise 

A(x) = 1, \x\ < a; =0, otherwise 

D(x) = H(Q), a < \x\ < 1 ; =0, otherwise (27) 

where a = sin Q = sin Q' and where 

H(Q) = if Q < tt/2, #(Q) = 2 if Q > tt/2 . 

The integral equations (p2]) and fl23|) become 

/oo 
K(x - a/)B(a/)(r(a/)<y, \x\ < a (28) 

-CO 

/OO /*OG 
iT(x - x')A(x')/(x')rfx + / K(x - x')D(x')t(x')dx' 
-OD J — OO 

/OO 
K 2 (x - x')B{x')a(x')dx' \x\ < B . (29) 

-oo 

Although these equations have to be solved in the stated intervals we can use 
their right sides to define their left sides for all real x. We define t(x) = for 
\x\ > 1. It is obvious that the extended equations have (unique) solutions 
if and only if the original ones do. Henceforth, we shall understand the 
functions / and cr to be defined for all real x. 
These equations read, in operator form, 

/ = t + KBa (30) 
a = KA f + KD t — K 2 Ba (31) 

where K is convolution with K and A, B, D are the multiplication operators 
corresponding to A, B, D (and which are also projections since A 2 = A, etc.). 

In view of the normalization requirements (|T9"D, the space of functions to 
be considered is, obviously, L l {\—a, a]) for / and L 1 ^— B, B}) for a. (L p is the 
p — th power integrable functions and L°° is the bounded functions.) Since 
K(x) is in L X (IR) fl L°°(R), it is a simple consequence of Young's inequality 
that the four integrals in (|28|) and ( ^9l) are automatically in L X (]R) fl L°°(R) 
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when / G L x ([— a, a]) and a G L 1 ^— B, B]). Particular, the integrals are 
in L 2 {R), which allows us to define the operators in ([3(]), ([H]) as bounded 
operators on L 2 (R). In addition, t is in but not in L 2 (R). To repeat, 

solutions in which / and a are in L 1 (M) automatically satisfy / — t and a 
are in L P (R) for all 1 < p < oo. 

THEOREM 1 (Uniqueness). The solutions f(x) and a(x) are unique 
and positive for all real x. 

REMARK: The uniqueness implies that / and a are even functions of x 
(because the pair f(—x),a(—x) is also a solution). The theorem implies 
(from the definition (ffi5|)) that cr(A) > for all real A and it implies that 
p(k) > for all \k\ < tt/2. It does not imply that p(k), defined by the 
right side of (22), is non-negative for all \k\ > tt/2. We shall prove this 
positivity, however, in Lemma 3. Note that the positivity of p is equivalent 
to the statement that f(x) < 2t(x) for all \x\ < 1 because, from (p2|) and the 
evenness or a, p(7r — k) = (1/tt) — p(k). 

Proof: By substituting (|30D into (|3T|) and rearranging slightly we obtain 



(1 + K 2 )a = K(A + D)t + K\\ - B)a + KAKBo . (32) 

Since K 2 is positive definite, 1 + K 2 has an inverse 1/(1 + K 2 ), which we can 
apply to both sides of ([32]). The convolution operator 

R = K(l + K 2 )- 1 (33) 

has a Fourier transform |sech(c<j/4). The inverse Fourier transform is pro- 
portional to sech(27rx) (see (|55|)), which is positive. In other words, R is not 
only a positive operator, it also has a positive integral kernel. 
We can rewrite (|32| ) as 

(l- W\<T = R(A + D)t = £ (34) 

with 

W = RK(1 -B) + RAKB = r\k~(1- A)KB\ . (35) 

The middle expression shows that the integral kernel of W is positive. 

Clearly, £ > as a function and £ G L 1 (IR)nL 2 (R). Also, W has a positive 
integral kernel. We note that ||i?|| = 1/2 on L 2 (IR) since y/(l + y) < 1/2 for 
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y > 0. Also, \\K\\ = 1, and ||1 - B\\ = 1, ||A|| = 1, ||S|| = 1. In fact, it is 
easy to check that ||iL4iifS|| < 1/2. From this we conclude that \\W\\ < 1 
on L 2 (M) and thus 1 — W has an inverse (as a map from L 2 (R) — > L 2 (M)). 
Therefore, we can solve ([32|) by iteration: 

o= (\ + W + W 2 + W Z + ...^ji . (36) 

This is a strongly convergent series in L 2 (M) and hence (|36| ) solves (^) in 
L 2 (IR). It is the unique solution because the homogeneous equation (1 — 
W)4> = has no solution. Moreover, since each term is a positive function, 
we conclude that a is a positive function as well. QED. 

Lemma 1 (Monotonicity in B). When B increases with Q fixed, cr(x) 
decreases pointwise for all x el. 



Proof: Since 1 — A is fixed and positive, we see from the right side of (35) 
that the integral kernel of W is monotone decreasing in B. The lemma then 
follows from the representation (|36|) . QED. 

Lemma 2 (Monotonicity in B). When B increases with Q fixed, f(x) 
increases pointwise for all x G K. This implies, in particular, that p(k) 
increases for all \k\ < it/ 2 and decreases for all n/2 < \k\ < n. 

Proof: Consider equation (^) for the case A = 0. Theorem 1 and Lemma 
1 hold in this case, of course. We also note that their proofs do not depend 
on any particular fact about the function Dt, other than the fact that it is 
a non-negative function. From these observations we learn that the solution 
to the equation 

(1 + K 2 B)S = Kg (37) 

has the property, for all x G M, that S(x) > and that S(x) is a non- 
increasing function of B, provided only that g{x) > for all iGl 

Another way to say this is that the integral kernel of V — (1 + K 2 B)~ 1 K 
is positive and is a pointwise monotone decreasing function of B. 



Now let us rewrite d37|) as 

K 2 (1 + K 2 )~ l g = Ug + K(l + K 2 )~ l (l - B) S (38) 

with 

U = KBV = KB(1 + K 2 B)~ l K . (39) 
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The operator U has a positive integral kernel since V, K, and B have one. 
As B increases the second term on the right side of ( j38|) decreases pointwise 
(because (1 — B) decreases as a kernel and S decreases, as we have just 
proved). The left side of fl38|) is independent of B and, therefore, the first 
term on the right side of (^) must increase pointwise. Since this holds for 
arbitrary positive g, we conclude that the integral kernel of U, in contrast to 
that of V, is a pointwise increasing function of B. 

Having established the monotonicity property of U let us return to /, 
which we can write (from ([3(])) as 



./ 



(l + UD)t + UAf 

[1 + UA + {UAf + ■ • • ] (1 + UD)t 



(40) 
(41) 



The series in ([41]) is strongly convergent (since ||y4|| = 1 and \\U\\ < 1/2) and 
thus defines the solution to (|4"0"D. Since U is monotone in B, ( flip tells us that 
/ is also pointwise monotone, as claimed. 



Equation fl26|) tells us that p(k) is increasing in B for \k\ < tt/2 and is 
decreasing in B for 7r/2 < \k\ < tt . QED. 

THEOREM 2 (Monotonicity in B). When B increases with Q fixed, 
N/N a and M/N increase. When B = oo, we have 2M = N, and when 
B < oo we have 2M < N (for all Q ). 

Proof: The integral for N/N a , in fli~9|) can be written as J^JAp + 
( l/27r) .D], and this is monotone increasing in B since p is monotone for 
\k\ < n/2 and A(k) = for \k\ > tt/2. 

If we integrate (^) from A = — oo to oo, and use the fact that J K = 1 
from (E4f), we obtain 



N 



/Q poo pB 

p{k)dk= / a(A)dA+ / a(A)dA (42) 
Q J-oo J-B 



p(k) dk 

which becomes, after making use of the normalization (|i~9|) 
l = 2f + § / +/ o-(A) dA 





poo 


J — oo 


Jb 



(43) 



Now the integrals in (^) decrease as B increases by Lemma 1 and converge 
to as B — > oo, while N/N a increases, as we have just proved. It follows 
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that M/N increases monotonically with B, reaching M/N = 1/2 at B = oo. 
If B < oo then M/N < 1/2 since a is a strictly positive function. QED. 

We turn now to the dependence of a, p on Q, with fixed B. First, Lemma 
3 (which was promised in the remark after Theorem 1) is needed. 

Lemma 3 (Positivity of p). For all B < oo, allQ < tt, and all \k\ < tt, 

we have p(k) > 0. 

Proof: As mentioned in the Remark after Theorem 1, the positivity of p is 
equivalent to the statement that f(x) < 2 t(x) for all \x\ < 1. We shall prove 
f(x)<2 t(x) here. 

Owing to the monotonicity in B of / (Lemma 2) it suffices to prove the 
lemma for B = oo, which we assume now. We see from (|4l|) that for any 
given value of a the worst case is Q > tt/2, whence H(Q) = 2 and D > 0. 
We assume this also. 

For the purpose of this proof (only) we denote the dependence of f(x) on 
a by f a (x). 

We first consider the case a = 0, corresponding to Q = tt. Let us borrow 
some information from the next section, where we actually solve the equations 
for B = oo, Q = tt and discover (Lemma 5) that f(x) < 2t(x) for \x\ < 1 
(for U>0). 

We see from ( |40|) or (|4T|) that f a is continuous in a and differentiable 
in a for < a < 1 (indeed, it is real analytic). Also, since the kernel 
K(x, y) = K(x — y) is smooth in (x, y) and t(x) is smooth in x G (—1, 1), it 
is easy to see that f a is smooth, too, for x G (—1, 1). Equation (BH) defines 
f a (x) pointwise for all x and / a (x) is jointly continuous in a,x. 

In detail, (|40|) reads 



Take the derivative with respect to a and set h a (x) = df a (x)/da. Observe 
that U does not depend on a. We obtain 




f a {x) = t(x) + 2 



J —1 J a J —a 



/ +/ U{x,x')t(x')dx' + U(x,x')f a (x')dx' . 



(44) 



h a {x) = U (x, a) [f a {a) - 2t{a)] + U (x, -a) [f a (-a) - 2t(-a 



)] 




(45) 
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(This equation makes sense because f a (x) is jointly continuous in x,a and 
t(x) is continuous for \x\ < 1. Recall that / and t are even functions of x. 
Note that U here is the kernel of with B = oo, i.e., U = K 2 (l + K 2 )-\ 
which is self-adjoint and positive as an operator and positive as a kernel.) 
Equation (25) can be iterated in the same manner as fl2T|) (since \\U\\ = 



1/2) 

h a (x) = [U + UAU + UAUAU + ■ ■ -}(x,a)F(a) = T(x,a)F(a) , (46) 

where [ ■ ](x,a) denotes the integral kernel of T = [ ■ ], and where F(a) = 
fa(a) ~ 2t(a) is a number. As an operator. T is self-adjoint and positive. 

Now U has a positive kernel and thus T(x, a) > 0, so h a (x) < for all x 
if and only if F(a) < 0. We have already noted that F(0) < 0. 

We can integrate (E^) to obtain 



f a (x) = f (x) + / h a i(x)da! = f (x) + / T(x, a')[f a >{x) - 2t(a')} da' . 
Jo Jo 

(47) 

If we subtract 2t(a) from this and set x = a, we obtain 

F(a) = G(a)+ [ T(a, a')F(a') da' , (48) 
Jo 

where G(a) = f (a)-2t(a) < 0. Another way to state (||) is F = G + TAF. 

Equation ( ^8|) implies that F(a) < for all a, as desired. There are 
two ways to see this. One way is to note that T is monotone increasing 
in a (as an operator and as a kernel), so T < U + U 2 + • • • = K 2 < 1. 
Therefore, (0) can be iterated as F = [1 + TA + TATA + •••]£?, and 
this is negative. The second way is to note that f a {a) (and hence F(a)) is 
continuous in a. Let a* be the smallest a for which F(a) = 0. Then, from 
= F(a*) = G(a*) + / Q a * T(a*, a')F(a')da' < 0, which is a contradiction. 

From F(a) < we can deduce that f a (x) — 2t(x) < for all \x\ < 



1. Simply subtract t(x) from both sides of (^7|). Then / a (x) — 2t(x) = 
{/o( a; ) — 2t(x)} + (TAF)(x). The first term { } < by Lemma 5, which we 
prove in Sect. 6 below, and the second term is < (since F < 0) QED. 

Lemma 4 (Monotonicity in Q). Consider the dependence of the solu- 
tion to ffity), ( (ffj[ ) on the parameter < a < 1 for fixed B < oo. For Q < n/2 
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(i.e., H(Q) = D = 0), both f and a increase pointwise as a increases. For 
Q > 7r/2 (i.e., H(Q) = 2,Dt = 2(1 — A)t), both f and a decrease pointwise 
as a increases. 

If, instead of the dependence on a, we consider the dependence on < 
Q < n of p{k) (which is defined by for all \k\ < it) and of a (A) (which 
is defined by (||3p for all real A), then, as Q increases, 

p{k) increases for < \k\ < n/2 and decreases for n/2 < \k\ < n 
cr(A) increases for all real A . (49) 

Proof: Concerning the monotonicities stated in the second part of the lemma, 
(|49l), we note that as Q goes from to tt/2, a increases from to 1, but when 
Q goes from tt/2 to 0, a decreases from 1 to 0. Moreover, H(Q) = in the 
first case and H(Q) = 2 in the second case. This observation shows that the 
first part of the lemma implies the statement about a in ([4*9]) . The statement 
about p in (|^) also follows, if we take note of the cos k factor in fl22|) . 

We now turn to the first part of the lemma. The easy case is Q < tt/2 or 
H(Q) = 0. Then ([|l]) does not have the UDt term and, since U has a positive 
kernel and since A has a kernel that increases with a, we see immediately 



that / increases with a. Likewise, from (p4|), (|35|) , we see that W and £ 
increase with a and, from (pop, we see that a increases. 

For Q > 7r/2 or H(Q) = 2, we proceed as in the proof of Lemma 3 by 
defining h a (x) = df a (x)/da and proceeding to ([R]) (but with U given by 
(39|)). This time we know that F(a) < (by Lemma 3) and hence h a (x) < 0, 
as claimed. The monotonicity of a(x) follows by differentiating ( |2"9"| ) with 
respect to a. Then (da(x)/da) = (VAh a )(x) + V(x, a)F(a), where V(x, y) is 
the kernel of V, which is positive, as noted in the proof of Lemma 2. QED 

THEOREM 3 (Monotonicity in Q). When Q increases with fixed 
B, N/N a and M/N a increase. When Q = ir, N/N a = 1 (for all B ), while 
N/N a <lifQ<n. 

Proof: From fl42|) , N/N a = 2 f^ B o~ + 2 f^o~ and this increases with Q 
by (H). Also, by (gD, N/N a = J Q Q p. When Q = tt, we see from (||) that 
J Q Q p = J^(1/2tt) = 1, so N/N a = 1. To show that N/N a < 1 when Q < n 
we use the monotonicicity of a with respect to B (Lemma 2) and Q (Lemma 
3) (with cr(A) = the value of cx(A) for B = oo, Q = n) to conclude that 
N/N a < 2 f B B cr + 2 J^°cr = J^ff = 1 — 2 'a < 1, since a is a strictly 
positive function. 
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Finally, from (22) we have that M/N a = J^ B <J, and this increases with 
QbyfH). QED 



6 Solution for the half-filled band 



In the case of a half-filled band, we have iV = N a ,M = M' = N/2 and, 
from Theorems 2 and 3, Q = tt, B = oo. In this case the integral equations 
([22] ) and (|23|) can be solved. We use the notation po(k) and cr (A) for these 
solutions. 



Substituting (p2|) into (|23|) where, as explained earlier, we use po(k) = 
1/27T in the first term on the right side of (|23|) . Then the integral equation (|23|) 
involves only 0"o(A) and can be solved by Fourier transform. Using equations 
(|24| ) it is straightforward to obtain the solution for a and its Fourier transform 
cr as 



ao(A)rfA 



2cosh([/a;/4) 



(50) 



/A . 1 f°° Jq{uj) cos (a; A) 

or (A) = — / : , .;, du. 



27r J cosh(a;[//4) 



where 

2 Z^ 2 1 P 

J (c<j) = — cos(u cos 9) d9 — — I cos(u sm9)d9 

Wo ft Jo 

is the zeroth order Bessel function. 



Next we substitute ([51]) into (|22|) and this leads (with (p4|) ) to 



(51) 



(52) 



,, s 1 cos A; f°° cos(a; sin k) Jn(^) , 
Po(k) = — + / K - ; T ,r ' du. 



2tt 7i 



1 + e" u / 2 



(53) 



The substitution of flSBI ) into Q2"5| ) finally yields the ground state energy, 
E , of the half-filled band as 



Er 



Na Na 

2 ' 2 



-4iV„ 



J (q;) Ji(cl>) 
u;(l + e^V 2 ) 



do; 



(54) 



22 



where Ji(u) = n 1 Jq sm(u ship) smpdp = ljtc 1 J* cos(a;sinp) cos 2 p dp is 
the Bessel function of order one. 

REMARKS: (A.) When there is no interaction (U = 0), K is a 5-function; 
we can evaluate (|5lD and (153) as 



°"o( A ) = - — ,] . q ? |A| < 1; =0, otherwise , 

Itx\J\ — A^ 

1 7T 

Po(&) = — ) 1^1 < tt! =0, otherwise . 

71 2 

This formula for po(fc) agrees with what is expected for an ideal Fermi gas. 

(B.) The U — > oo limit is peculiar. From (pUD we see that a (0) = 1/2, 
so J a = 1/2, but from ( |51~1) we see that 0" O (A) — ► in this limit, uniformly 
in A. On the other hand po(k) — > 1/2tt, for all \k\ < tt, which is what one 
would expect on the basis of the fact that this 'hard core' gas becomes, in 
effect, a one-component ideal Fermi gas of N = N a particles. 

We now derive alternative, more revealing expressions for a, p. 

For (Tq (A) we substitute the integral representation (|52"D for J into (|5TD 
and recall the Fourier cosine transform (for a > 0) 

r^x) du= rp j i (55) 

Jo cosh(u;a:) \2a>y cosh(7rx/2a;) 

Then, using 2 cos a cos 6 = cos(a — b) + cos(a + 6) we obtain 



de 

> o . 



cosh[2vr(A + cos #)/£/] 



(56) 



An alternate integral representation can be derived similarly for po{k), but 
the derivation and the result is more complicated. We substitute (l + e^) -1 = 
J2^=i(~ 1)™ exp[— nx], with x = uU/2, into (^) and make use of the identity 
2], 6.611.1] (with a = -is ± c in the notation of |Za 6.611.1]) 



roc 

2 / e-™ Jq(uj) cos(cus) du = [(-c - is) 2 + l] ~ 1/2 + [(c - is) 2 + l] ~ 1/2 
Jo 

(57) 
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for c > 0, and where the square roots [ 
positive real part. This leads to 



-1/2 



111 



are taken to have a 



1 cos k 

"" (k) = + 



oo 

J2(-l) n+1 1 [(-nIJ/2 - i sin kf + l] 



•1/2 



71=1 



+ [{nU/2 -isink) 2 + l] 



-1/2 



(58) 



We can rewrite the sum of the two terms in ( pq) as a single sum from 
n = — oo to oo, after making a correction for the n = term (which equals 
cosfc/| cos/c| for fc 7^ 7r /2). We obtain 

cos /c 



1 


1 + 


cos k 




cos/c 


2^ 


1 


1 + 


cos k 




cos/c 


27 



2tt 

n= 

cos k 1 



^ (-l) n [{nU/2-isinkf + 1] 

=— oo 



-1/2 



7T 



isin/c) 2 + 1 sm(7T^J 
(59) 



The contour C encompasses the real axis, i.e., it runs to the right just below 
the real axis and to the left just above the real axis. 

The integrand has two branch points y± on the imaginary axis, where 
V± — {2i/U)(sink ± 1). In order to have the correct sign of the square root 
in the integrand we define the branch cuts of the square root to extend along 
the imaginary axis from y + to +oo and from y_ to — oo . We then deform 
the upper half of the contour C into a contour that runs along both sides 
of the upper branch cut and in two quarter circles of large radius down to 
the real axis. In a similar fashion we deform the lower half of C along the 
lower cut. As the radius of the quarter circles goes to oo this gives rise to 
the following expression. 



Po(k) 



1 

2^ 



1 + 



cos k 
I cos/c I 



cos k 
2ttU 



[/_(*) + !+(*)] >° 



(60) 



where 



I ± (k) 



da 



i±sinfc sinh(27ra/[/)A/(a: =F sin k) 2 — 1 



(61) 
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By introducing the variable a = cosh x ± sin k we finally obtain the simple 
expression 



sinh {(27r/[/)(coshx ± sin k)} 

As a consequence of expressions (|60|) and (^) for p, we have the crucial 
bound needed as input in Lemma 3: 

Lemma 5 (p bounds). When B = oo , Q = it, and U > 

1/2tt < po(fc) < V 71 " for 0<|A;|<7r/2 

< p (£:) < 1/2tt for 7r/2<|A;|<7r (63) 

Equivalently , fo(x) < 2t(x) for all \x\ < 1. 

Proof: When 7r/2 < | k \ < n and cos k < the first term [ ] in (|BT)D is 
zero while the second term is positive (since I±(k) > 0). On the other hand, 
when < \k\ < n/2, Theorem 1 shows that po(k) > 0. Thus, we conclude 
that po(k) > for all \k\ < it. 

From (|22]) and the positivity of a we conclude that po(k) < 1/2-k when 
7r/2 < \k\ < n. From the positivity of po(k) when n/2 < \k\ < n we conclude 
that the integral in (^) is less than 1/27T for all values of < sin k < 1 and, 
therefore, l/2vr < p (k) < l/n for < \k\ < tt/2. QED 




7 Absence of a Mott transition 

A system of itinerant electrons exhibits a Mott transition if it undergoes a 
conducting-insulating transition when an interaction parameter is varied. In 
the Hubbard model one inquires whether a Mott transition occurs at some 
critical U c > 0. Here we show that there exists no Mott transition in the ID 
Hubbard model for all U > 0. 

Our strategy is to compute the chemical potential p + (resp. p_) for 
adding (resp. removing) one electron. The system is conducting if p + = p_ 
and insulating if p + > p_. 

In the thermodynamic limit we can define /i by p = dE(N)/dN, where 
E(N) denotes the ground state energy with M = M' = N/2. As we already 
remarked, this choice gives the ground state energy for all U, at least in the 
thermodynamic limit. 
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The thermodynamic limit is given by the solution of the integral equa- 
tions, which we analyzed in Sect. |5|. In this limit one cannot distin- 
guish the odd and even cases (i.e., M = M' = N/2 if N is even or 
M = M' - 1 = (N - l)/2 if N is odd.) and one simply has M/N = 1/2 in 
the limit iV a — > oo. In this case Theorem 2 says that we must have B = oo. 
Then only Q is variable and Theorem 3 says that Q is uniquely determined 
by N provided N < N a . 

In the thermodynamic limit we know, by general arguments, that E(N) 
has the form E(N) = N a e(N/N a ) and e is a convex function of N/N a . It 
is contained in (E5|) when N/N a < 1. A convex function has right and left 
derivatives at every point and, therefore, /i + = right derivative and //_ = left 
derivative are well defined. Convexity implies that /i_ < 

For less than a half-filled band it is clear that fi + = /x_ since E(M, M) is 
smooth in M = iV/2 for iV < N a . The chemical potential cannot make any 
jumps in this region. But, for N > N a we have to use hole-particle symmetry 
to calculate E(N). The derivatives of E(N), namely fi + and can now be 
different above and below the half-filling point N = N a and this gives rise to 
the possibility of having an insulator. We learn from (§) that 

/i+ + //_ = U, (64) 

and hence /i + > /i_ if /i_ < C//2. 

We calculate //_ in two ways, and arrive at the same conclusion 



oo 



^ 1 ^ 7 cu[l + exp(^/2)] 



(65) 



The first way is to calculate /i_ from the integral equations by doing per- 
turbation theory at the half-filling point analyzed in Sect. ^|. This is a 



'thermodynamic' or 'macroscopic' definition of \X- and it is given in Sect. |7T 
below. (From now on /i_ means the value at the half-filling point.) 

In Sect. [F72] we calculate /x_ 'microscopically' by analyzing the Bethe 
Ansatz directly with N = N a — 4 electrons. Not surprisingly, we find the 
same value of This was the method we originally employed to arrive at 
§ Eq. 23]. 

Before proceeding to the derivations of (|65|), we first show that (65) im- 



26 



plies yU_ < U/2 for every U > 0. To see that < U/2 we observe that 



Mo) 
Mo) 







a; 



Ji(u) du = -. 



(66) 
(67) 
) is in 



Then < U/2 holds if fi'!_(U) < 0, which we turn to next. Here, 
H 6.561.17] and fl67|) is in ||, 6.511.1] 

Expanding the denominator in the integrand of ([88]) and integrating term 
by term, we obtain 



/i_(E0 = 
using which one obtains 
t>"-(U) 



2D- 1 : 



2 

2vri 



1 + 



~2~J 



n 



n 2[/2N3/2 



-v2 



7T 



C 



[/2 z 2x3/2 



Sin 7TZ 



dz 



(68) 



where we have again replaced the summation by a contour integral with the 
contour C encompassing the real axis. The integrand in (^) is analytic 
except at the poles on the real axis and along two branch cuts on the imagi- 
nary axis. This allows us to deform the path to coincide the imaginary axis, 
thereby picking up contributions from the cuts. This yields 



If 



(ya-l)^ sinh(^) 



dy < , for all U > . (69) 



Thus, we have established n+{U) > fi~(U), and hence the ID Hubbard model 
is insulating for all U > 0. There is no conducting-insulating transition in 
the ground state of the ID Hubbard model (except at U = 0). 



7.1 Chemical potential from the integral equations 

As we noted, we take B = 00 and Q < tt. In fact we take Q = 71 — a with 
a small. (In the notation of Sect. [5], a = sinQ, but to leading order in a, 
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sinQ = 7i — Q and we need not distinguish the two numbers.) Our goal is 
to calculate 5E, the change in E using (|25|) and oW, the change in N using 
(|T9|); p_ is the quotient of the two numbers. 

As before, we use the notation p{k) for the density at Q = tx — a and 
Po(k) for the density at Q = n, as given in (|53[) , (|60|). 

We start with N. As explained earlier, p — 1/2tx is odd around 7r/2 so, 
from (|19D, 

/•« z-Tr-a 

p = 2 p = 2 p + 2 — 



N a J-q Jo Jo Ir, 2-k 

p + -(tt - 2a) w 1 + 2a f p (0) - - 1 . (70) 



7T \ 7T 



In the last expression we used the fact (and will use it again) that p is contin- 
uous in k and a (as we see from ([41])); therefore, we can replace J Q a p by ap (0) 
to leading order in a. We learn from (|70|) that 5N/N a = 2a (po(0) — < 0. 
The calculation of SE is harder. From (p5|) 

— = — 4 / p cos k = —4 / p cos k — A I p cos 

A a Jo JO J a 

pit— a 2 rn—a 

~ — 4apo(0) — 4 / (p )cosA; / cosfc 

J a 27T 7T ,/ a 

/■tt/2 1 

= -4ap (0)-8 / (p-— )cosk 

Ja 27T 

tt/2 



= — 4ap (0) H (1 — sin a)— 8 / pcosfc + 8 / pcos/c 

2tt Jo Jo 

4a Z"^ 2 4 / >7r ^ 2 
« +4ap (0) 8/ 5pcosfcH 8/ p cos/c, (71) 

77 Jo ^ Jo 

where 5p = p — po- The last two terms in ([71]) are the energy of the half- filled 
band, N = N a . 

Our next task is to compute 5p to leading order in a. It is more conve- 
nient to deal with the function 8f = / — fo and to note (from (|26|) ) that 
Jq^ 2 p(k) cos kdk = Jq 1 f(x)\/\ — x 2 dx. We turn to ( f4T[) and find, to leading 
order, that 

/ « (1 + 2C/) t - f/lt + 2£M£?t = /o + C/A/o - 2£Mt (72) 
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with f Q = (1+2U) t. We note that U = K 2 (l + K 2 )~ l since B = oo (see ((J|)) 
and has a kernel which we will call u(x —y). If g is continuous near (in our 
case g = f or g = t) then (UAg)(x) = f" a u(x - y)g{y) dy « 2au{x - 0)0(0) 
to leading order in a. We also note from (|26|) that /o(0) = Po(0). Therefore, 



tt/2 



<5pcos k m a 



Po(0) - - 



7T 



vl — x 2 u(x) dx . 



(73) 



The integral in (|73|) is most easily evaluated using Fourier transforms and 
Plancherel's theorem, 



VT^e luJX dx 



and from 



cos(tt;x) Vl — <ix 



2 / cos(cj sin 5) cos 2 6d6 = - J x (u) (74) 
Jo u 



u\x)e 



x dx = [1 



jU/2\ 



By combining these transforms we can evaluate 8E from (|7T|). 



5E 

— =2a 

N„ 



p (0) - - 

7T 



2-4 



w[l + exp(o;C//2)] 



By dividing ( |76|) by ( |70D we obtain (|65|) . 



(75) 



(76) 



7.2 Chemical potential from the Bethe Ansatz 

The evaluation of the chemical potentials from the Bethe Ansatz is reminis- 
cent of the calculation of the excitation spectrum of the ID delta-function 
Bose gas solved by one of us |RJ . We consider the case of a half-filled band. 
To use our results in the previous sections, which hold for M, M' odd, we 
calculate /i_ by removing 4 electrons, 2 with spin up and 2 with spin down, 
from a half-filled band. This induces the changes 



iV -> iV-4 



M -> M - 2 = N a /2 - 2 (77) 
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The equations (15) and (|T6|) determining the new k' and A' now read 

M-l 

Nak'j = 2tt/; + 0(2 sin k] - 2A' p ) , j = 3, 4, • • • , N - 2 (78) 

AT-2 M-l 

^ #(2sin fcj - 2A'J = 2rcJ' a - ^ 6(A' a - AJ,) , a = 2, • • • , M - 1, 

3=3 (3=2 

(79) 

Under the changes (|7?D, the values of /' and J' are the same as those of / 
and J, namely, 

I', = Ij, J' = 3,4,... ,N-2 
J' a = J a , a = 2, 3, • • • M - 1 

so they are centered around the origin with k' total = fctotai- 

The removal of four electrons causes the values of k and A to shift by 
small amounts, and we write 

kf j = k j + w w(k j ), A' a = A a + — u(A 8 ) 

By taking the differences of ( fTSf) and ( |l5|) , and (^) and (|T6|) , and keeping 
the leading terms, one obtains 



u>(%) = — ^ d> ( 2 sin k J ~ 2A /3) [ 2 cos k M k j) ~ 2u ( A ?)] > ( 80 ) 



/3=2 
iV-2 



4#(2A a ) - — Y 0'(2A a - 2 sin kj) [2u(A«) - 2 cos kjw(t 



3=3 
M-l 



a (3=2 

In deriving these equations we have used facts from our analysis of the 
integral equations, namely that when M = M', — Ai = Am ~ oo (i.e., = oo 
in the limit N a — > oo) and that when N = N a , —k\ = k^ ~ — &2 = fcjv-i ~ ^ 
as AT a — ■> oo. Without using these facts there would be extra terms in (p0|) 
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and ([SID , e.g., 6 (2 sink j — 2Ai) + 9(2sinkj — 2A M ), which is « because 
-Ai = Am w oo. 

By replacing the sums by integrals and making use of P0|) and (|2T|), we 
are led to the coupled integral equations 

/oo 
^(sinA;- A)s(A) dA (82) 
-OO 

/7T 
K(sin/c — A)r(k) cos kdk 
-IT 

/OO 
K 2 (A - A')s(A')dA' (83) 
-oo 

where 

r(k) = w(k) Po (k) , S (A) = u(A)<7 (A) . (84) 



Equations (82) and (^) can be solved as follows. Note that the third 
term on the left side of (|83|) vanishes identically after substituting (p^ ) for 
r(/c). Next introduce the Fourier transforms (E3) and 



oo 



e 4wA 0(2A) dA = - f — ) e-^l^ 4 , (85) 



and we obtain from fl83| ) 

/IS 2 f°° sinwA , . . 

7r Jo cosh(u;cy/4j 



Thus, from (|82| 



4 r sin(o;sinfc) ^ 
W ~W w (l + ^) ' (87) 



Note that we have r(—k) = —r(k) and s(— A) = — s(A). 
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The chemical potential /x_ for a half-filled band is now computed to be 
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i=3 



r(/c) sin k dk 
o w(l + e^/ 2 ) 



da; 



(88) 



which agrees with 



8 Conclusion and acknowledgments 

We have presented details of the analysis of ground state properties of the 
one-dimensional Hubbard model previously reported in |J. Particularly, the 
analyses of the integral equations and on the absence of a Mott transition 
are new which have not heretofore appeared in print. 

It is important to note that in order to establish that our solution is indeed 
the true ground state of the ID Hubbard model, it is necessary to establish 



the existence of ordered real solutions to the Bethe Ansatz equations (14) 
and, assuming the solution exists, proofs of (a.) and (b.) as listed at the end 
of Sect. 3. The fulfillment of these steps remains as an open problem. 

We are indebted to Daniel Mattis for encouraging us to investigate the 
jump in the chemical potential as an indicator of the insulator-conductor 
transition. We also thank Helen Au-Yang and Jacques Perk for helpful dis- 
cussions. FYW would like to thank Dung-Hai Lee for the hospitality at the 
University of California at Berkeley and Ting-Kuo Lee for the hospitality 
at the National Center for Theoretical Sciences, Taiwan, where part of this 
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